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Abstract. We construct the first X-ray selected sam- 
ple of broad line radio-loud AGN from the EMSS sur- 
vey. In order to test unifying schemes for radio-loud ob- 
jects, their spectral and statistical properties (both flat 
and steep spectrum objects) are examined and compared 
with those of other samples of blazars. The X-ray selec- 
tion allows us to explore properties of radio-loud quasars 
10-100 weaker in the radio band than classical samples. 
The most convincing interpretation of our results is that 
there are no radio-loud quasars whose synchrotron emis- 
sion peak reaches the EUV-soft X-ray band at these 
(radio) flux levels. Moreover, due to the comparatively 
weak non-thermal emission, a quasi-thermal component 
appears to contribute at optical-UV energies. The detec- 
tion of sources at low radio fluxes also reveals a large 
population of steep spectrum quasars and the lack of the 
predicted turnover in the quasar (radio) counts. The evo- 
lution of X-ray selected radio-loud quasars does not sig- 
nificantly differ from that of radio-selected ones, and flat 
spectrum and steep spectrum sources appear to behave 
quite similarly. 

Key words: Galaxies: active - Galaxies: evolution - 
Galaxies: luminosity function - Quasars: general - Radio 
continuum: galaxies - X-rays: galaxies 



1. Introduction 

The selection energy band in the construction of complete 
samples constitutes one of the main biases and a source of 
confusion in the understanding of the Active Galactic Nu- 
clei (AGN) phenomenology, which has often led to useful 
but somehow unphysical classifications. 

Among radio-loud objects, a specific and instructive 
example - which will constitute a reference in the follow- 
ing - is given by the case of BL Lac objects. Samples of 
BL Lacs selected in the radio and X-ray bands revealed 
systematic differences in their spectral energy distribu- 
tions (SED), polarization, variability properties and evo- 
lution (e.g. Urry & Padovani 1995). Different models have 
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been proposed to account for this apparent "dichotomy" 
(Stocke et al. 1985; Ghisellini & Maraschi 1989; Celotti et 
al. 1993), and only recently it has been explicitly recog- 
nized that an intrinsic continuous distribution in the po- 
sition of the peak of the broad band SED of these sources 
might naturally lead to the detection of different types 
of objects depending on the selection band (Padovani & 
Giommi 1995). This new perspective has allowed us to 
have a better view of the intrinsic, physical properties of 
BL Lacs and speculate on their origin. 

The SED of all blazars can be well represented by two 
main spectral components, believed to be produced as syn- 
chrotron and inverse Compton emission, respectively. One 
suggested possibility then is that the broad band spec- 
tral distributions of 'X-ray selected' (XBL) and 'radio- 
selected' (RBL) objects are simply distinguished by the 
energy location of the peak of the synchrotron emission: at 
UV/X-ray in the former ones (High peak BL Lacs, HBL) 
and around the optical in the latter ones (Low peak BL 
Lacs, LBL; Padovani & Giommi 1995), while the Comp- 
ton peak is located in the 7-ray band (from MeV to TeV 
energies). 

Furthermore, the entire blazar population, ranging 
from HBL, to LBL to Flat Spectrum Radio Quasars 
(FSRQ) can be interpreted, so far, as a sequence from line- 
less, low power sources, to broad lined, powerful quasars 
characterized by the (decreasing) peak energies of the two 
spectral components, the increasing source radio power 
and the increasing ratio between the Compton and syn- 
chrotron luminosities (Fossati et al. 1998). This trend in 
the SED can be physically accounted for by an increase 
in an external soft photon field along the sequence, in 
agreement with the observed behavior (Ghisellini et al. 
1998). Clearly, this scenario makes definite predictions on 
the relationship between (radio) source power and spectral 
shape. 

In order to understand the properties of a population 
it appears therefore crucial to select sources in different 
spectral bands, and compare their characteristics. To this 
aim, i.e. disentangling the intrinsic features from the selec- 
tion induced ones, we have considered radio-loud quasars 
- which so far have been extracted from radio surveys - 
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and studied their broad band and cosmological properties 
when selected in X-rays. 

The sample comprises the radio-loud quasars (i.e. 
broad line AGN) detected in the Einstein Medium Sensi- 
tivity Survey (EMSS; Gioia et al. 1990, Stocke et al. 1991). 
More specifically, we have constructed the first statisti- 
cally complete X-ray selected sample of (39) radio-loud 
quasars, in order to: a) determine whether their spectral 
properties differ from those of the radio-selected quasars 
and other classes of radio-loud objects; b) discuss the 
properties of these objects in the context of the beaming 
scenario and emission mechanisms for blazars; c) deter- 
mine the statistical properties of this new sample, in the 
light of the unification scheme for powerful radio-loud ob- 
jects (e.g. Urry & Padovani 1995 and references therein), 
where flat spectrum radio quasars are high power (FR II) 
radio galaxies in which the relativistically beamed compo- 
nent dominates. 

With the selection criteria adopted, the sample com- 
prises both Flat Spectrum (FS, blazar-hke) and Steep 
Spectrum (SS) radio quasars [|. It should be reminded 
that the radio emission in the latter ones is believed to 
originate from the extended radio components. In these 
sense the SS quasars will be not directly included in the 
comparison of blazar properties, but we will study and 
examine their optical. X-ray spectral characteristics and 
their statistical properties with respect to the FS quasars 
in this sample. 

After discussing the sample selection (Sect. 2), we 
present the broad band spectral properties of the objects, 
distinguishing between flat and steep radio sources (Sect. 

3) and compare them with those of other blazars (Sect. 

4) . The statistical and evolutionary properties of the new 
sample are derived in Sect. 5, while in the final Sect. 6 we 
discuss our results. 

Hereafter we adopt cosmological parameters Hq — 50 
km s^^ Mpc~^ and qq — 0. 

2. The sample 

The EMSS provides us with well defined X-ray selected 
samples of extragalactic sources. The virtually complete 
identification status makes it a very powerful tool to study 
statistical properties of classes and sub-classes of astro- 
nomical sources. 

Radio information for all the (extragalactic) EMSS 
objects is available from VLA snapshots at 6 cm taken 
during the identification process (e.g. Maccacaro et al. 
1994), with a flux limit of about 1 mJy, corresponding 
to 5 times the rms of the observations. Since we intend to 
compute the radio spectrum using the recently released 
data at 20cm from the NRAO VLA Sky Survey (NVSS; 

^ Note that the sample is not selected on the point-like opti- 
cal appearance, and thus the SS sources possibly include Broad 
Line Radio Galaxies. 



Condon et al. 1998), we cut the EMSS sample at ^ > -40°. 
With this cut, that excludes only 5% of all the AGNs and 
only one radio-loud, the EMSS contains 79 radio detected 
AGN, of which 43 have a two-point spectral indexQ be- 
tween radio and optical (aro) larger than 0.35 (Delia Geca 
et al. 1994) and are therefore classifled as radio-loud. We 
exclude another four objects for the reasons detailed in 
the Appendix, thus constructing a complete sample of 39 
radio-loud quasars, that we dub the EMSS Radio-Loud 
quasar sample (ERL). 

In the following we will refer to FS and SS spectrum 
radio-loud objects, where the division is formally consid- 
ered to be at ar=0.7 (see e.g. Perlman et al. 1998) which 
leads to a similar statistics for the two sub-classes (20 
FS and 19 SS). As the distribution of a^- does not show 
any suggestion of being the sum of two different (bimodal) 
distributions (see Fig. |l|) , the choice of a dividing value is 
somewhat arbitrary. In any case, the results of this work 
are not affected by a different choice (the case ar=0.5 has 
been also explored). 

The list of sources and information is presented in Ta- 
ble 1^. Columns are as follows: (1) Name of the source in 
the EMSS catalog - numbers refer to notes to the Ta- 
ble, indicating either other names from literature or the 
source radio morphology; (2) Redshift; (3) 5 GHz flux from 
EMSS, unless otherwise stated, in mJy; (4) 1.4 GHz flux 
from NVSS; (5) Radio spectral index; (6) Magnitude, from 
EMSS; (7) unabsorbed X-ray flux (0.3-3.5 keV band), from 
EMSS; (8) Galactic iVn, from Stark et al. (1992); (9) X- 
ray spectral index from EMSS data; (10) X-ray spectral 
index from ROSAT data (Brinkmann et al. 1997); (11-13) 
and ttrx- 

The distributions of radio spectral indices and lumi- 
nosities in the three bands are plotted in Fig. |] for FS 
and SS, respectively. 

3. Spectral properties 
3.1. X-ray spectral indices 

The X-ray spectral indices for the sample have been esti- 
mated from the IPC hardness ratios (see e.g. Maccacaro 
et al. 1988), assuming a power law spectral shape and 
Galactic line of sight absorbing columns, and compared 
with those determined by ROSAT. The former ones span 
a large interval —2.2 < ax < 2.4, and within the un- 
certainties significantly correlate with the ROSAT values 
(see Fig. ||). Unfortunately the large errors do not allow 
either to detect any difference between FS and SS sources, 

^ While the original EMSS two point spectral indices Ofro, 
"ox, Orx were estimated at 5 GHz, 2500 A and 2 keV, comoving 
energies, here we recompute them at 5 GHz, 5500 A and 1 keV, 
comoving energies, in order to compare with values for other 
samples. For the K-correction in the optical and X-ray bands 
it is assumed Uo — = 1 (see e.g. Delia Ceca et al. 1994), 
where the spectral indices a are defined by F^, ~ 
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Fig. 1. Distributions of radio spectral indices and lumi- 
nosities for the ERL. The FS and SS subsamples are indi- 
cated by the dotted and continuous line, respectively. The 
axis on top of the figure refers to the spectral index ar- 




Fig. 2. The ROSAT spectral index (from Brinkmann ct 
al. 1997) vs. the EMSS spectral index derived from the 
Hardness Ratio (from Maccacaro et al. 1988); FS and SS 
sources are indicated as crosses and circles, respectively. 



or to infer the origin of the X-ray component (e.g. steep 
synchrotron dominated versus flat Compton dominated 
emission) . 

No relationship has been found between the X-ray and 
radio spectral indices, the source fluxes or redshifts. 



3.2. Broad band properties 

In order to quantify the total power and SED properties of 
the sources and compare them with those of other samples, 
we compute the luminosities in the radio, optical and X- 
ray bands (Lr, io, L^) and the corresponding broad band 
spectral indices (aro, ctox, Q^rx)- 

Significant linear correlations are found between the 
luminosities in the three bands. However, since a possible 
culprit for this kind of correlation is the common redshift 
dependence, we test its influence by using a partial corre- 
lation analysis (e.g. Kendall & Stuart 1979). The partial 
correlation coefficients are defined as: 



ri2_3 = 



(ru - ri3 X r23) 



13) X i^-rh) 



(1) 



where ri2, ria, and are the Spearman rank test coef- 
ficients for two variables at a time. By assigning to index 
"3" the redshift variable and interchanging L^, and 
Lo as variables "1" and "2" we obtain the following re- 
sults: r„_2 = 0.441; rox.z = 0.325; rro.z = 0.381. The nuU 
hypothesis is rejected at the 95% level for the last two oc- 
currences, while the radio and X-ray luminosity seem to 
be correlated at the 99% level. 

No differences (according to the Kolmogorov-Smirnoff 
test) are found between the distributions relative to FS 
and SS, neither in luminosities (probability p > 20%; see 
Fig. ^ nor in broad band spectral indices (p > 47%; see 
Fig. 0) nor in the trends, suggesting a behavior of the com- 
pact beamed component independent of the large scale 
radio flux dominance. 

The shape of the SED appears however to be clearly 
related to the radio luminosity: in fact, the only statisti- 
cally signiflcant trends found are of both and arx with 
Lr (see Fig. Although the correlation might be due to 
the common dependence of the two spectral indices on 
itself, its relevance is supported by two facts: a) L,. does 
not correlate with aox; b) L,. spans the largest interval in 
luminosity (about flve decades compared to the three cov- 
ered by Lo and L^). In any case, the presence of such a 
trend is indicative of the fact that through X-ray selection 
we detect objects with a rather limited range in optical lu- 
minosity but a large span in the radio one, extending the 
range of sampled Lr towards faint sources compared to 
radio-selected quasars (see Sect. ^). 

On the other hand, it is interesting to consider the 
relative behavior of the broad band spectral indices, as 
shown in Fig. ^. arx, which spans the smallest interval, 
and aox, which spans the largest one, do not show a strong 
correlation (98% level, according to the Student's-t test), 
while (a part for a conspicuous object, the weak lined 
MS0815. 7-1-5233 - see also Sect. 4) aro decreases for in- 
creasing aox (at > 99% level). Interestingly, this trend is 
qualitatively consistent with a change in the frequency po- 
sition of the peak of the synchrotron spectral component. 
Furthermore, by comparing Figs. || and ^, it is clear that 
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Table 1. The sample 

Name z F| F^^^ mag Qrosat a, 



MS0012 


5-0024 


1 


695 


24.0 


14.9 


-0 


40 


18 


51 


2 


41 


2 


9 


1 


cn+O.TO 
^"-0.63 











32 


1 


33 


0.67 


MS0038 


8-0159^ 


1 


675 


322.0 


1034.6 





97 


16 


86 


3 


07 


2 


9 


-0 


qi +0.69 
^'--1.17 











53 


1 


54 


0.88 


MS0039 


2-0206 


1 


083 


5.6 


7.0 





19 


19 


00 


1 


48 


2 


9 




— 











30 


1 


33 


0.66 


MS0053 


3-1035 





308 


27.6 


40.6 





32 


18 


44 


3 


33 


3 


7 















42 


1 


29 


0.72 


MS0136 


3+0606 





450 


16.3 


67.5 


1 


18 


18 


52 


2 


72 


3 


7 


2 


-0.78 
c:n+0.50 
^^^-0.63 


1 


78±0 


16 





41 


1 


30 


0.72 


MS0152 


4+0424 


1 


132 


184.1 


323.8 





47 


19 


17 


6 


28 


4 


2 














63 


1 


07 


0.78 


MS0225 


5-1052 


1 


038 


60.0 


132.2 





66 


19 


34 


1 


56 


2 


5 


-0 




1 


08+0 


15 





56 


1 


29 


0.81 


MS0226 


8-1041^ 





620 


96.2 


293.2 





93 


18 


32 


2 


53 


2 


5 


1 


^ O+0.47 
^■^-0.48 





92±0 


19 





54 


1 


35 


0.82 


MS0232 


5-0414^ 


1 


450 


522.9 


1495.9 





87 


16 


28 


7 


19 


2 


5 


1 


1 ir+0.38 
^"^-0.39 


1 


18±0 


06 





51 


1 


49 


0.85 


MS0311 


8-0801 


1 


250 


8.2 


70.3 


1 


78 


19 


12 


3 


38 


5 


4 


1 


in+0.92 
^'J-1.19 











44 


1 


17 


0.70 


MS0402 


0-3613* 


1 


417 


2019.0 


1132.3 


-0 


48 


17 


17 


33 


98 





8 










75±0 


16 





59 


1 


12 


0.78 


MS0438 


6-1050 





868 


26.8 


51.2 





54 


18 


97 


3 


48 


5 


9 


-1 


45+1^90 











46 


1 


19 


0.71 


MS0521 


7+7918 





503 


22.1 


102.6 


1 


28 


17 


45 


5 


52 


7 


9 





qo+1.41 











35 


1 


33 


0.69 


MS0808 


0+4840 





700 


75.0 


160.3 





63 


17 


96 


2 


97 


4 


6 


1 


4C+0.80 


1 


32±0 


16 





48 


1 


36 


0.78 


MS0815 


7+5233^ 





624 


1.5 


3.6 





73 


20 


75 


3 


15 


4 


3 


2 


,-,0+0.93 
'-'•-'-0.72 











37 





94 


0.56 


MS0822 


0+0309 





577 


31.5 


85.3 





83 


17 


94 


1 


46 


3 


5 


2 


oc+9.99 
°^^-1.78 











41 


1 


49 


0.78 


MS0828 


7+6601 





329 


98.8 


235.3 





72 


20 


21 


4 


90 


4 


2 


-0 


U«-0.71 











69 





95 


0.78 


MS0833 


3+6523** 


1 


112 


369.0'' 


1275.7 


1 


03 


18 


16 


3 


84 


4 


4 


-0 







46±0 


24 





65 


1 


29 


0.87 


MS0849 


0+2845 


1 


273 


336.0 


232.6 


-0 


31 


19 


09 


1 


73 


3 


2 




— 











62 


1 


29 


0.85 


MS0850 


2+2825 





922 


37.0 


68.0 





51 


18 


27 


1 


37 


3 


2 


-2 


48+2^33 











42 


1 


46 


0.78 


MS0958 


4+6913'^ 





928 


23.0 


59.8 





79 


19 


08 


2 


68 


3 


7 


-0 


9/1+0.67 
^^-1.04 





24±0 


18 





47 


1 


22 


0.73 


MS1003 


6+1300^'* 





648 


125.0 


305.7 





74 


17 


45 


2 


75 


3 


8 


1 


00+0.68 
JO_0.70 











49 


1 


46 


0.82 


MS1050 


9+5418 





995 


167.1 


266.0 





39 


19 


25 


4 


89 





8 





^9+0.49 
^^-0.59 


-0 


02±0 


45 





63 


1 


12 


0.79 


MS1138 


6+6553^ 





805 


132.8 


506.4 


1 


11 


18 


30 


1 


26 


1 








-7C-+0.53 
'•-'-0.56 


1 


00±0 


27 





57 


1 


48 


0.89 


MS1234 


9+6651^ 





860 


142.2 


335.2 





71 


17 


99 


5 


39 


1 


9 


-1 


2^ + 1.71 

— OO 











53 


1 


29 


0.79 


MS1311 


1+3210^ 





303 


7.0 


132.0 


2 


44 


18 


86 


3 


29 


1 


1 


-0 


q-i +1.09 











39 


1 


24 


0.69 


MS1326 


6+2546^'^ 





986 


45.6 


89.3 





56 


17 


80 


2 


40 


1 


1 





7f.+0.47 
'D-0.50 





60±1 


00 





41 


1 


45 


0.77 


MS1340 


7+2859 





905 


250.0 


252.1 





01 


17 


05 


4 


88 


1 


2 





r r + 0.98 
"^•-'-1.41 


1 


06+0 


38 





47 


1 


44 


0.80 


MS1442 


8+6344 


1 


380 


445.0 


689.6 





36 


17 


19 


2 


93 


1 


7 


1 


r)c+0.47 
'^•-'-0.47 


1 


12+0 


16 





53 


1 


51 


0.87 


MS1623 


4+2712^° 





525 


240.0* 


531.0 





66 


18 


41 


7 


78 


3 


4 


1 


qo+0.32 
•^°-0.32 


1 


96+0 


39 





62 


1 


16 


0.80 


MS 1640 


0+3940 





540 


31.0 


54.4 





47 


18 


25 


4 


39 


1 








Qq+0.41 
^•^-0.39 


1 


03±0 


05 





41 


1 


29 


0.72 


MS1657 


1+3524 





949 


5.6 


5.2 


-0 


06 


20 


24 


1 


11 


1 


7 







1 


30±0 


11 





39 


1 


21 


0.67 


MS1837 


8+4538^-^1 





958 


210.0 


801.4 


1 


11 


19 


38 


3 


60 


5 


9 





84+.79 











70 


1 


12 


0.84 


MS2037 


3-0035^ 





609 


23.1 


73.1 





96 


18 


39 


2 


21 


6 


6 


1 


Qq+0.68 
^•^-0.67 
09+1.25 
^^^-3.03 


1 


30±0 


43 





42 


1 


34 


0.74 


MS2113 


8+0455 


1 


001 


122.0 


316.8 





79 


20 


37 


2 


53 


6 


5 


-0 





39±0 


00 





71 


1 


02 


0.82 


MS2141 


2+1730^^ 





213 


417.0 


651.8 





37 


15 


85 


8 


94 


8 





2 




1 


44±0 


23 





46 


1 


49 


0.81 


MS2143 


2+1429 


1 


387 


50.3 


226.2 


1 


25 


20 


40 


1 


91 


7 


















67 


1 


06 


0.80 


MS2247 


8-0703 


1 


710 


44.7 


82.4 





51 


20 


82 


3 


82 


3 


6 


1 


9f.+0.59 
^"-0.62 











63 





90 


0.72 


MS2329 


3-3827^'^ 


1 


195 


670.0 


543.6 


-0 


17 


16 


20 


3 


08 


1 


5 


-1 


05+2^02 











46 


1 


65 


0.87 



Flux @ 5 GHz in mjy; Flux @ 1.4 GHz in mjy; " Unabsorbed flux in (0.3-3.5 keV) in xlO"*^ erg cm^^ s^^ Hydrogen 
column density in xlO^°cm"^; " from Gregory & Condon (1991); ' from Becker et al. (1991). Notes: (1) 4C02.04 (2) Triple 
radio source (3) PKS 0232-04, double radio source (4) PKS 0402-362 (5) Weak-lined, see text (6) 3C 204 (7) Radio 
core+extended (8) WE 1138+65 (9) WAT (10) B2.2 1623+27A (11) OU+462 (12) OX 169 (13) PKS 2329-384, HB 



this spectral behavior is connected with the source radio 
luminosity, as the flatter is Oro, the lower is L^. Notice also 
that the trend does not simply relate to the optical and 
X-ray luminosities, which therefore appear to be worse 
tracers of systematic behaviors in the SED shape. 

We conclude that: a) the ERL sample (and the FS and 
SS subsamples) significantly extends the radio luminosity 

range (but not the optical one) of radio-loud quasars; b) 
the broad band spectral indices behavior is consistent with 



being caused by a change in the synchrotron peak energy; 
c) the trend in the SED shape correlates with the radio 
luminosity. 

4. Comparison with other samples 

We now compare the spectral properties of the ERL sam- 
ple with those of radio-selected quasars (FSRQ and Steep 

Spectrum radio-loud quasars, SSQ) and (radio and X- 
ray selected) BL Lac objects and examine our findings in 
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Fig. 3. Distribution of arx, Qox, ctro as a function of Lr- 
SS are represented by circles and FS by crosses. The linear 
regression for the entire sample is plotted for arx and a^o 
(the separate regressions for FS and SS are consistent with 
the total), while vs is not statistically significant. 



the light of the blazar spectral sequence scenario already 
described and the unification of powerful radio sources. 
Does the X-ray selection allow us to discover powerful 
lined sources whose synchrotron emission peaks at high 
(UV- X-rays) energies, analogously to X-ray selected BL 
Lac objects? Does the detection of ERL force us to re- 
consider the above scenario? Do the nuclear properties of 
steep spectrum objects differ from those of the flat spec- 
trum sources or their nuclear properties are independent 
of the extended emission? 

Let us first examine the former issue. 

The plane aox versus has been often used to char- 
acterize the SED of the above classes of sources. As men- 
tioned in the introduction, the position of the different 
objects in this diagram can be accounted for by the dif- 
ferent location of the synchrotron peak in their energy 
distribution. 

In Fig. ||a we show the ERL and other samples (see 
the figure caption) of both FSRQ and radio and X-ray 
selected BL Lacs. For clarity, we also plot in Fig ||b a 
schematic representation of the loci occupied by the dif- 
ferent samples, in which the mean and the la width of 
the distribution in the two variables are plotted. The ERL 
occupy the transition region (in the "boomerang" shaped 
blazar distribution) between FSRQ and X-ray selected 
HBL) BL Lac objects, roughly overlapping with radio- 
selected (~ LBL) ones. Note that the location of LBL and 
FSRQ in the spectral index plane is due to the dominance 
in the X-ray band of a flat Compton component. 



Fig. 4. Broad band spectral indices: arx (upper panel) 
and aro (lower panel) vs. aox for the ERL. Symbols like 
in Fig. I 

Within the spread of the above quantities, this indi- 
cates that even X-ray selection is unable to detect radio- 
loud broad-lined objects with high peaked synchrotron 
components. Intriguingly, the object in the left-bottom 
corner, in the "HBL region", is the already mentioned 
MS08I5.7-h5233, a weak fined AGN, BL Lac-fike (see 
Stocke et al. 1991, Tables 11 & 13 and notes therein). 

Furthermore, the broad band spectral indices arx and 
ttro of the ERL tend to be lower (Aa ^0.1-0.2) than 
those of the FSRQ, with the fiattest sources coinciding 
with the very same objects which have lower Lr- In other 
words ERL resemble radio selected quasars in the optical 
and X-ray bands, while are comparatively weaker in the 
radio one, indeed suggesting a larger contribution (peak) 
to the SED at higher energies. This result is fully consis- 
tent with the high energy/X-ray selection of the sample. 

It is also interesting to explore the relationship among 
the luminosities of sources belonging to different samples, 
as the lack of ERL with SED resembling those of HBL 
might be due, within the above scenario, to their relatively 
higher luminosity. In Fig. ^ the above samples are plotted 
in the Lx vs Lr plane (luminosities uncorrected for evolu- 
tion). While in the X-ray band the ERL are on average 
a factor of ^ 10 weaker than the FSRQ and similarly the 
SSQ, in the radio (core) luminosities they are of the order 
of 100 (10) times weaker than the FSRQ (SSQ) (e.g. Fos- 
sati et al. 1998; Maraschi & Rovetti 1994). The comparison 
also shows that the Lx of ERL covers the same range spun 
by HBL, although it extends at even higher Lx- Neverthe- 
less there is no indication of a trend for the weakest ERL 
to have SEDs more resembling highly peaked sources, as 
their radio luminosity range is above that of HBL. We 
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stress that in the very same EMSS sample HBL are in- 
deed found in large number (Morris et al. 1991, Wolter et 
al. 1991), excluding observational biases against detection 
of ERL peaking in the X-ray band. 

These findings are therefore in global agreement with 
the expectation of the blazars sequence scenario described 
in the introduction, thus re-enforcing the view of a strong 
connection between the spectral energy distribution and 
the radio luminosity and suggesting again that the latter 
might be a good indicator of the source power. 

However, it should be noted that ERL and LBL span 
a similar range of three decades in L^. While one would 
expect some overlapping along the sequence - as also sug- 
gested by the existence of sources changing from being 
BL Lac-like to quasar-like and vice versa, and RBL with 
FR II extended morphology - the extent of the common 
luminosity range seems to suggest more of a parallel be- 
havior, rather than a continuous trend, of lined and line- 
less sources. However, it appears that at lower the 
ERL might be filling the intermediate region between LBL 
and HBL. These two facts, and thus the overall distri- 
bution of radio and X-ray luminosities of ERL, might 
be be ascribed to the increasing dominance of a quasi- 
thermal optically-UV component in quasars of increas- 
ingly lower (non-thermal) radio power. Indeed the typi- 
cal a-co and aox of the low radio power ERL (see Fig. ||) 
corresponds to SED whose peak is located in the optical- 
UV region. The presence of a quasi-thermal component 
would indeed flatten a-co and steepen ctix in sources with 
decreasing non-thermal continuum. This 'excess' contri- 
bution could be identified with both/either 'blue bump' 
emission in quasars (FS) and/or contamination from the 
host galaxy in SS (possibly Broad Line radio galaxies). 



5. Statistical and evolutionary properties 

Let us now examine the statistical properties of the sam- 
ple within the framework of the unification scheme for 
powerful radio sources. More specifically, do the steep ra- 
dio sources evolve consistently with the postulate of the 
unification scheme (i.e. similarly to fiat sources)? 

5.1. Number counts 

The number counts of radio-loud AGN of the EMSS have 
been presented by Delia Ceca et al. (1994). Here we con- 
sider any possible distinction related to the radio spectral 
index (FS vs SS), and compare the ERL counts to those 
of other samples of radio-loud objects. 

The (integral) logA^ ~ log 5 distributions for FS, SS, 
all ERL and BL Lacs are shown in Fig. |^, and the para- 
metric results derived from their linear fits are reported 
in Table ^ (The fit to the counts is done by using the 
individual points in a differential form, see e.g. Gioia et 
al. 1991 for a description of the procedure). 
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Fig. 5. a) Broad band spectral index a^o vs. aox for the 
ERL (filled triangles and circles for FS and SS, respec- 
tively) compared to: EMSS BL Lacs (empty circles). Slew 
BL Lacs (empty triangles), IJy BL Lacs (empty stars), 
2Jy FSRQ (empty squares). Just for graphical purposes, 
we have not included in the plot the most recent blazars 
selected by Perlman et al. (1998) and Laurent-Muehleisen 
et al. (1999). b) The lower panel shows the loci of the dif- 
ferent samples, the black circle indicates the ERL sample, 
while the triangles the BL Lac and FSRQ samples. The 
errorbars represent the width (Icr) of the distributions in 
ttro and ttox respectively. 



The first point to be noticed is the marginal (though it 
extends over the whole range in flux) evidence for flatter 
counts of FS, which dominate at the highest flux levels. 
The count distributions of FS is in complete agreement 
with what found for all AGN - mostly radio quiet- in 
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Table 2. Statistical properties 





ERL 


FS 


SS 


Count slope " 


1.9 [1.7-2.2] 


1.8 [1.6-2.1] 


2.0 [1.7-2.3] 




0.92 ± 0.40 


0.99 ± 0.40 


0.85 ± 0.39 




0.78 ± 0.05 


0.76 ± 0.06 


0.81 ± 0.07 


c 




6.6 [5.6-7.4] 


6.4 [5.6-7.1] 


7 




2.9 [2.5-3.1] 


3.0 [2.7-3.2] 


LF" 


1.6 [1.4-1.7] 


1.6 [1.3-1.8] 


1.9 [1.6-2.2] 


logK^ 


-8.74 


-8.95 


-9.00 



In square brackets the 1 a confidence range 
Power-law index of a linear fit (integral slope) 
± lo" values 

if is in units of yipc~^ L^l~^ , where L44 is the luminosity in units of 10'*'' erg s 
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Fig. 6. Monochromatic (@1 keV) X-ray vs monochro- 
matic (@5 GHz) radio luminosity of the ERL, EMSS BL 
Lacs, Slew BL Lacs, IJy BL Lacs, 2Jy FSRQ (symbols as 
in Fig. H). The dotted line indicates the "nominal separa- 
tion" between HBL and LBL in the luminosity plane (as 
defined in Fossati et al. 1998). 



comprises 34 SSQ and 50 FSRQ 0. The (poorly defined) 
count distribution in the radio band is consistent with be- 
ing Euclidean and also with those of FS and SS (at 90%) 
and of ERL (at 95%) in the X-ray band. 

The most interesting aspect is the extension to a fac- 
tor 100 lower fluxes provided by the X-ray selection of the 
ERL. Indeed assuming typical arx (~ 0.8), (ar) and (ax), 
an X-ray flux level of ~ 10~^^erg cm~^ s^^ corresponds 
to ~ 0.3 Jy in the radio band (both for FS and SS). At 
these levels, the counts of FS exceed the prediction of the 
beaming model proposed by PU92. Furthermore if one ex- 
trapolates to lower radio fluxes (say by a factor ~ 10 as in 
the X-ray band), the radio counts of X-ray selected ob- 
jects are marginally consistent with the Euclidean extrap- 
olation from higher fluxes and thus largely exceed (factor 
>10 for FS and ~3 for SS) those predicted by the beaming 
model, which drop below ~ 0.1 Jy. Note that because of 
the X-ray threshold the radio counts of ERL can be only 
a lower limit to the actual surface density of radio-loud 
quasars. 

It is also relevant to compare the ERL counts with 
those of BL Lacs extracted from the EMSS. As shown 
in Fig. 7, the BL Lac distribution is flatter, dominat- 
ing by a factor ~ 10 at 10^^^ erg cm^^s^^, and an 
'inversion' of the two populations occurs just above ~3 
X 10"^'^ erg cm~^ s~^. Therefore while at high X-ray fluxes 
RL quasars are as rare, or rarer, than BL Lacs, deep, small 
area surveys are bound to find large numbers of them. 



the EMSS (1.61 ± 0.06, Delia Ceca et al. 1992), while the 
slopes of SS and all ERL are marginally steeper (but only 
at 1 cr level) than those of all AGN. Clearly the statistics 
does not allow us to perform more complex fits, although 
a single power law might be inadequate to represent the 
distributions. 

How do the ERL counts compare with those of FSRQ 
and SSQ? Padovani & Urry (1992, hereafter PU92) ex- 
amined the statistical properties of the best defined radio 
sample, namely the 2 Jy (Wall & Peacock 1985). This 



5.2. Redshift distribution 

All sources are below redshift of ~ 2, with a broad peak 
around z ~ 1. This is due to the sky coverage of the 
EMSS, combined with the LF of the AGN. There is only 

^ We note here that the average radio spectral indices of the 
two subsamples of FS and SS are fully consistent with those 
of FSRQ and SSQ once the same definition is adopted (more 
precisely we obtain: FS: 0.32 ± 0.36 (0.15± 0.34), SS: 0.96 ± 
0.14 (0.85 ± 0.20), for a dividing = 0.7(0.5) 
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Fig. 7. Integral number counts for FS (crosses), SS (cir- 
cles) and all ERL (filled circles). The lines represent the 
corresponding linear fits: heavy line for ERL, dashed line 
for FS, dot-dashed line for SS. The flatter distribution 
(squares) represents the EMSS BL Lacs (from Wolter et 
al. 1991). Results of the fits are reported in Table 0. 




Fig. 8. Redshift distribution for SS sources (solid line) and 
FS sources (dotted line). 



marginal evidence for SS to be at redshifts lower than FS 
(See Table ||. Further, the probability of getting by chance 
a larger KS value is 19.6%). 

The typical z are globally similar (within 1 a) to the 
average redshifts of FSRQ and SSQ selected in the radio 
band (2Jy), while they are significantly larger than those 
of FR II radio galaxies, as already noted by PU92. 

5.3. Fe/K 

In order to quantify the evolutionary properties of the 
sample, we estimated, again subdividing it into FS and 
SS, the Vc/Va parameter (Avni & Bachall, 1980), a gen- 



eralization of the V/V„iax statistics (Schmidt 1968) when 
complete samples at different flux limits are combined. 
The results are shown in Table ^ Positive evolution is 
required for the subsamples and the whole population at 
> 4(7 level, while no statistical difference can be found 
between the FS and SS behavior. 

A systematically lower Vc/Vn (though consistent 
within ~ 2a) is found for the quasars of the 2 Jy sample 
(PU92), possibly indicating a tendency of faster evolution 
of the X-ray selected objects. ^ 

In order to compare the results with those relative to 
the 2Jy objects, we characterize the evolution in terms 
of an exponential pure luminosity function of the form 



Liz) 



-Ct 



where 1/C is the look-back time in 



units of the Hubble time. Note however that for the whole 
AGN sample in the EMSS, this form was found to be 
adequate only if a luminosity dependent value for C was 
introduced (Delia Ceca et al. 1992). The derived evolution 
rate is consistent (at ~ 2a, see Table ||) with that inferred 
for FSRQ (and once again it is indistinguishable between 
FS and SS). As our sample is limited to z < 2 nothing can 
be deduced on a possible decline in the number density of 
sources at high redshifts (e.g. Dunlop & Peacock 1990). 

Also a power law luminosity evolution L(z) = 
L(^2^Q-^{l+z)^ can well represent the cosmological changes. 
The required dependence on z is similar to what is found 
for radio-quiet objects (Table ^). The results for the FS 
sources are again consistent with the whole EMSS AGN 
behavior, while SS require a slightly higher value of 7 (at 
la). 

For computing fluxes and luminosities we have as- 
sumed that — 1 - supported by the fact that {a^) ~ 1 
in the ERL sample and by a study of the IPC Hard- 
ness Ratio for the EMSS AGN (Maccacaro et al. 1988). 
Note however that for X-ray spectra with a slope differ- 
ent from 1 the effect of the cosmological evolution has 
the same functional form of the K-correction factor and, 
in first approximation, the best-fit value of 7 is lower (if 
ttx > 1) or higher (if < 1) by an amount ~ |1 — ax 
(for |1 - Oxl < 1)- 

5.4. Luminosity function 

Finally, we built the luminosity functions (LF) of FS and 
SS, both as observed and applying the two evolutionary 
models to de-evolve luminosities at z = (Figure ||a, b). 
As expected, no difference is found between flat and steep 
spectrum objects (the extension of FS a factor of ~ 10 
more in luminosity is due to only one source). The distri- 
butions have been fitted with power-laws (see results in 
Table ^) , as the large uncertainties do not allow to better 
characterize their shape. 



^ As this behavior is opposite to that of BL Lacs derived 
from the EMSS, this indirectly (re)confirms the reality of the 
uncommon negative evolution found for XBL. 
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Fig. 9. Differential luminosity functions for the SS and 
FS sources, in step of 0.5 logL (symbols as in previous 
figures; filled circles and stars indicate de-evolved SS and 
FS sources respectively) . The SS points have been slightly 
shifted artificially on the luminosity axis for clarity. The 
de-evolving form is: upper panel exponential, lower panel 
power law. 



At odds with these results, a significant difference in 
the number density of flat and steep spectrum objects is 
found in radio selected samples (2 Jy), with FSRQ numer- 
ically dominating above ~ 10'''^ erg s~^. Converted into 
the X-ray band, this luminosity corresponds to ^ 10^^ 
erg s~^, where however a smaller number of SS, within 
the volume explored, are found. 

The slope of the local LF of ERL is significantly steeper 
than that of FSRQ, revealing a first clear difference be- 
tween the population of radio selected and X-ray se- 
lected flat spectrum quasars. This discrepancy is of course 
present also with respect to the X-ray LF computed by 
Maraschi & Rovetti (1994) from the radio LF of FSRQ 
and SSQ over the common range in X-ray luminosity. 

The total number densities of FSRQ is predicted to be 
- 2.2 Gpc-3 above - 5 x lO*^ gj.g (PU92) with 
a parent population of 81 FR II Gpc~'^. Locally we find 
twice this density. Also, SS are ^ 7 times those predicted 



by extrapolating below the cutoff the PU92's model pre- 
dictions: in other words we do not find evidence for either 
a cut off or a flattening of the LF below X-ray powers of 
~ 3 X lO''^^ erg s^^. XBL have a slope slightly flatter than 
the ERL one and a local number density (at say 10'*'' erg 
s~^) about 50 times larger. 

Also the ratio of FS/SS found in the EMSS signifi- 
cantly differ from that estimated by Maraschi & Rovetti 
(1994), i.e. a ratio of 30 between the two populations, while 

- at the corresponding X-ray luminosity - we find similar 
numbers of FS and SS. In particular the density of FS is 
about one tenth than their predictions. ^ 

The same 'excess' of SS with respect to FS in the X- 
ray selected sample clearly reflects in the relative ratio 
of SS/FS=19/20 found in the EMSS, compared to 34/50 
found in radio selected survey at higher flux levels (2 Jy) 

- despite SS do not dominate at low radio luminosities. 
This seem to imply that the radio selection is more biased 
than the X-ray one in finding FSRQ with respect to SSQ. 

5.5. Predictions: number counts and redshift distributions 

Finally, let us consider the number counts and redshift 
distributions predicted by the above findings. In Fig. 10a 
the number counts extrapolated at flux levels observable 
by Chandra and Newton- XMM are reported. These have 
been derived by the integration of the evolving LF (of FS 
and SS) up to redshift 5 (between = lO'*^ — 

10'*'^ erg s"^). The hmits of the proposed Newton-XMM 
surveys are plotted in Fig. 10a. Given the areas that will 
be probably covered at each flux limit (^1,5, and 150 sq. 
deg. respectively in the faint, medium and bright survey), 
we expect about 2000 (300), 150 (50), 100 (45) SS (FS) 
quasars in each, if the measured evolution is taken at face 
value and the extrapolation holds at the faint fluxes. These 
samples will be large enough to allow refined prediction 
on the radio-loud quasars properties. The corresponding 
redshift distributions at different flux limits are shown in 
Fig. 10b, where SS can be seen to dominate at lower z at 
all flux thresholds. 



6. Summary and discussion 

We have constructed the first complete X-ray selected 
sample of radio-loud quasars. This has allowed us to in- 
vestigate both the spectral properties and possible trends 
among blazars and the statistics of such sources in the 
frame of unification models for radio-loud objects. 

Through X-ray selection (at the level of ^ 10"^"^ erg 
cm^* s^^) we explore the same range of redshifts and X- 
ray and optical luminosities of radio-selected objects, but 
we find quasars significantly fainter in the radio band. 



^ Note that for a flat spectral index Or ~ the threshold 
flux would be comparable at higher selection radio frequencies, 
therefore implying that the discrepancy is likely not due to the 
(radio) selection band. 
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Fig. 10. a) Predicted number counts for the FS (solid line) 
and SS (dashed line) at flux levels reachable by Chandra 
and Newton-XMM, in the soft band, for a population ex- 
tending up to Zmax 5. b) Rcdshift distributions of FS 
and SS, derived from the evolving LF, at different flux 
limits. 

As far as the spectral properties of blazars are con- 
cerned, FS ERL do not differ significantly from radio se- 
lected sources. However, the decrease in radio luminosity 
correlates with the increasing frequency of the spectral 
peak, as indicated by the trend of the broad band spec- 
tral indices. This supports the existence of a sequence in 
the blazar family, where the peak decreases in energy with 
increasing source radio power (Fossati et al. 1998; Ghis- 
ellini et al. 1998). 

Furthermore, although ERL are selected at X-ray lu- 
minosities comparable to those of HBL, we find no evi- 
dence for the existence of quasars whose synchrotron emis- 
sion component peaks at X-ray frequencies. This is also 
indicated by the similarity of of ERL and FSRQ, which 
suggests that the synchrotron component does not dom- 
inate the X-ray band luminosity (see e.g. Fossati et al. 
1998; see also Padovani 2000). 

On one hand in any unifying scenario some overlapping 
between objects of different characteristics might be nat- 



urally expected (likely corresponding to nuclei changing 
character or mismatching the extended radio morpholo- 
gies). On the other hand the ERL spectral distributions 
appear to strongly resemble those of LBL over a large 
range of powers (Lr spans three decades), where both 
broad lined and lineless objects co-exist. This suggests 
that the line typology might not be univocally related to 
the source power and broad band spectral shape. Indeed, 
one might envisage that the (occasional) appearance of 
broad line components in BL Lacs is related to a change 
occurring in the nuclear source - but at such low luminos- 
ity levels that the lines never become quasar-like - while 
the blazar sequence of power vs continuum spectral be- 
havior would be a separate effect and the resemblance of 
the broad band spectral indices of LBL and ERL might 
be indeed attributed to a similar trend of the synchrotron 
peak position with (radio) luminosity. 

However, from the paucity of information we have on 
ERL, this similarity can be due to the progressive domi- 
nance of a broad thermal component (blue bump or con- 
tribution from the host galaxy in FS and SS, respectively) 
over the non-thermal emission in objects of lower non- 
thermal radio luminosity. Accounting for this effect in 
fact reduces the intrinsic overlapping of the two popu- 
lations: decreasing the optical luminosity by considering 
only the (actual) non-thermal emission corresponds to 
move ERL towards FSRQ in the spectral index plane. If 
so, the tighter relation between source power and SED 
shape (and line strength) would be re-enforced by the 
ERL sample. Clearly this ambiguity can be solved by the 
estimate of the optical and X-ray spectra of ERL. Support 
to this view also comes from the findings of the REX sur- 
vey (Caccianiga et al. 2000), where a significant number of 
sources are found with clear optical excesses with respect 
to FSRQ and LBL (see their Fig. 6), while the radio and 
X-ray luminosities are simply linearly correlated. 

We stress here that also the strong difference in the 
statistical properties of lineless and lined samples indi- 
cates that the two populations are not simply analogous 
manifestations of the same phenomenon, where only the 
line intensities are different. 

We find no indication of any dependence of the spectral 
properties on redshift, nor any difference between the nu- 
clear properties of flat and steep spectrum radio quasars, 
when selected in the X-ray band, except possibly flatter 
counts for the FS. In other words the nuclear component 
seems not to be directly dependent on the extended (ra- 
dio) properties. 

From the statistical point of view, the radio and X- 
ray selection sample the same population at the brightest 
fluxes. The main new information we obtain is that at the 
lower flux levels sampled through the X-ray selection the 
radio source counts do not drop - and not even flatten - as 
predicted by some beaming models (e.g. PU92). Surveys 
at low radio flux levels mJy) will be able to establish 
the presence or absence of such a decline in the population. 
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The ERL sample has provided also with the possibil- 
ity of directly comparing the X-ray evolution of quasars 
with that of BL Lacs. This is particularly relevant both 
for understanding the reason of the negative evolution of 
BL Lacs in this band and in particular as a direct test 
of models postulating a physical evolution of radio-loud 
quasars into BL Lacs (e.g. Cavalicrc & Malquori 1999). 

The ERL population significantly evolves, at a rate in- 
distinguishable from that of FSRQ. With respect to the 
radio selection, the X ray one provides also with a differ- 
ent number ratio of flat and steep radio sources, with the 
latter comparable to the former ones. This number ratio 
does not appear to depend on the radio luminosity, but 
significantly depends on the radio - and possibly on the 
X-ray - flux, being the X-ray selection less biased at low 
radio fluxes. Both deeper radio surveys and the definition 
and study of optically selected blazars - as that which will 
be provided by the EIS surveys (e.g. Cagnoni & Celotti 
2001) - should confirm this trend. 

Due to the dramatic effect of beaming in altering the 
shape of the luminosity functions, it is not possible to 
simply use the above ratio as an indicator of the beaming 
angles subtended by FS and SS. Recent indications of the 
existence of a slower moving 'layer' jet component domi- 
nating the observed emission in radio galaxies (e.g. Laing 
1993, Chiaberge et al. 2000) further modify the expected 
number ratio of beamed and parent population sources. 
Furthermore as the luminosity functions do not show a 
turn over at low powers it is not even possible to provide 
model independent estimates on the total number of ob- 
jects of the two populations. Nevertheless it is clear that 
the beaming parameters and number density of the parent 
and beamed populations have to be revised in the light of 
the new and less biased information provided by the ERL, 
affecting the predictions on the radio-activity in AGN and 
galaxies as well as the contribution of these objects to the 
7-ray background. 

Appendix A: Sample selection: notes on 
individual objects 

Two of the 43 radio-loud objects of the EMSS have been 

classified as Narrow Emission Line objects (MS1252.4- 
0457 as a Seyfert 2 galaxy, and MS1414.8-1247 as an Hll 
region; Boyle et al. 1995) and therefore we remove them 
from the sample under study. 

MS0234.8-F0655 is detected at 6cm, but not at 20cm. 
By using the EMSS radio flux the source is radio- 
loud (aro=0.43), but with the NVSS flux limit (= 2.5 
mJy) ar < -0.81. Therefore both MS0234.8-h0655 and 
MS2134.0-F0028 (a^ = -0.88) have a highly inverted spec- 
trum and are removed from the sample. 

We examined also the case of MS1253. 6-0539 that has 
a very high upper limit from the EMSS data (5.5 mJy at 
5 GHz) and therefore in principle could be a radio source. 
However, the object is also not detected at 1.4 GHz by 



the NVSS at a limit of 2.5 mJy. Furthermore, it appears 
in Table 10 of Stocke et al. (1991) as having an ambiguous 
spectrum due to the poor signal to noise. We therefore do 
not include it in the sample. 

Note that MS0833.3-F6523 (a.k.a. 3C204) has F^^^ ~ 
30 mJy (from the EMSS VLA data) whfle Ftot = 369 
mJy: the extended flux dominates the emission and we use 
the latter to compute the radio spectral index. The same 
applies to MS1623.4-I-2712, for which the EMSS F^^™ = 43 
mJy while Ftot = 240 mJy. 
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